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The rotation of lo predicted by the Poincare-Hough model 

Benoit Noyelles^ 

NAmur Center for Complex SYStems (NAXYS), University of Namur, Dpt of 
Mathematics, 8 Rempart de la Vierge, B-5000 Namur, Belgium 

ABSTRACT 

This note tackles the problem of the rotation of lo with the 4-degrees of 
freedom Poincare-Hough model. lo is modeled as a 2-layer body, i.e. a triaxial 
fluid core and a rigid outer layer. 

We show that the longitudinal librations should have an amplitude of about 
30 arcseconds, independent of the composition of the core. We also estimate the 
tidal instability of the core, and show that should be slowly unstable. 

Subject headings: lo - Resonances, spin-orbit - Rotational dynamics - Interiors 
- Celestial mechanics 

1. Introduction 

Thanks to the Galileo space mission in the Jovian system, we have now interesting 
clues on lo's internal structure. In particular, its mass and its second-order gravity field 
coefficients J2 and C22 are known with a good accu racv (lAnderson et al.l 120011). Moreoy er, 



internal heating due to an inten se tidal dissipation (jPeale et al. 
expected to induce a fluid core ( ICassen et al.lll982l ). 



1979 



Lainey et al.ll2009l ) is 



We recently published a theo retical explora tion of the Poincare-Hough model applied 



to a synchronous 



dently by iHough 



rota ting b odv (iNoye 



( 18951) and iPoincarg ( 



led I2 0121) . This model, originally proposed indepen- 



19101 ) before being put in an Hamiltonian form by 



Touma fc WisdomI ( 120011 ) . describes the rotational dynamics of a triaxial body composed of 
a rigid mantle and an ellipsoidal cavity filled by an inviscid fiuid of constant uniform density 
and vorticity. This is a 4-degrees of freedom conservative model in which the core-mantle 
interactions result in pressure coupling at the core-mantle boundary. 



^email: Benoit.Noyelles@fundp.ac.be, also associated with IMCCE, CNRS UMR 8028, Paris Observatory, 
UPMC, USTL, 77 avenue Denfert-Rochereau, 75014 Paris, France 
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This paper proposes an application of this model to a realistic lo under the gravitational 
forcing of Jupiter. We build internal models and ap ply the Poincare-H ough model to them, 
and we use complete orbital ephemerides, here from (JLainey et al.ll2006l ). to model accurately 
the gravitational torque of Jupiter. After computation of lo's rotation we will discuss the 
problem of the elliptical instability of the fluid filling the core. 



2. The internal structure of lo 

We here consider lo to be an ellipsoidal body, its moments of inertia being < A < B < 
C. A is the moment of inertia with respect to lo's long equatorial axis, while C is related 
to the polar axis. The core is ellipsoidal as well, its principal axes of inertia being collinear 
to the ones of the whole body, the moments associated being written Ac, B^ and Cc- The 
known gravity and shape parameters of lo are gathered in TabO 



Table 1: Gravity and shape parameters of lo. / has not been measured independently of the 
others, it is deduced from the measured values of J2 and C22 combined with the hydrostatic 
equilibrium condition. 



Quantity 


Value 


Reference 


Mean density p 


3, 527.8 ±2.9 kg/m^ 


Anderson et al. ('2001) 


J2 


(1.8459 d= 0.0042) x lO'^ 


Anderson et al. ('2001) 


C22 


(5.537 ±0.012) X 10-^ 


Anderson et al. ('2001) 


1/ {MR'') 


0.37685 ±0.00035 


Anderson et al. ('2001) 


Mean radius R 


1,821.49 km 


Archinal et al. (2011) 


Subplanetary equatorial radius a 


1,829.4 km 


Archinal et al. (2011) 


Along orbit equatorial radius b 


1,819.4 km 


Archinal et al. (2011) 


Polar radius c 


1,815.7 km 


Archinal et al. (2011) 



The Poincare-Hough rotational model requires 5 internal structure parameters: 



Polar flattening ei 



2C-A-B 



2C 



J- 



MB? 



2 Q 1 



Equatorial ellipticity €2 



B-A 
2C 



2C. 



22- 



MB^ 



C ' 



■■^Oc ^c -Pc 



Polar flattening of the core €3 — ^„ 
Equatorial ellipticity of the core €4 



Be- Ac 
2Cc 



Relative inertia of the core 5 = ^ 
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From the definitions of J2 = (2C — A — B) / {2M R^) and tlie moment of inertia / = 

{A + B + C)/3, we liave 

C / 2 

^ = T + - J2 = 0.37808. (1) 

We can see from tliese data that the 2 parameters ei and 62 can be straightforwardly 
derived, yielding: 



ei 

£2 



4.88230 X 10^ 
2.92901 X 10 



-3 



The 3 parameters related to the size and shape of the core depend on its composition. 
We consider 2 end-members: either the core is made of pure iron, its density pr beiri g 
8, OOOfc^f/m^, or it is a eutectic mixture of FeS, yielding Pr = 5, 150kq/m^ (lUsselmanl Il975l ) . 
Another uncertainty is in the composition of the crust. lAnderson et al.l (J200ll ) considered 



2 possibilities: either a thin (< 50km) crust with a low density (< 2,600kg/m^), or a 
small amount of low density crust (that we propose to neglect) overlying a thicker (100- 
200 km) melt-rich asthenosphere (density between 3, 000 and 3, 200A:q/m^). T he reader can 



find additional information on lo's internal structure in iMoore et al 
assumptions, we get 6 interior models of lo (Tab|2]). 



torn . With these 



Table 2: Physical parameters of our 6 models. The models 1, 3 and 5 assume a pure iron 
core, while we have a eutectic FeS core in the models 2, 4 and 6. The models 3 and 4 also 
consider a 30 km-thick crust, while the models 5 and 6 have a 150 km-thick crust, ps is the 
density of the crust, and pm the one of the mantle. 



Models 


Pc 


Pm 


ps 


es 


£4 


6 




ikg/m^) 


{kg/m^) 


(kg/m^) 


xlOO 


xlOO 


XlOO 


1 


8,000 


3,291 


- 


1.722 


2.210 


1.654 


3 


8,000 


3,377 


2,400 


1.836 


2.673 


1.270 


5 


8,000 


3,409 


3,100 


1.877 


2.819 


1.188 


2 


5,150 


3,243 


- 


0.980 


1.057 


6.503 


4 


5,150 


3,337 


2,400 


1.029 


1.243 


5.175 


6 


5,150 


3,351 


3,100 


1.033 


1.253 


5.160 



We can see that the 6 models can be splitted into 2 groups, the discrimination coming 
fr om the compos i tion o f the core. The equations used to derive these models can be found 



m 



Novelles et al.l (120 llh . Eq.l to 12. 
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3. Rotational dynamics 

As already said, we use the Poincare-Hough model to represent the rotational dynamics 
of lo. This model considers a rigid outer layer, composed of the mantle and eventually the 
crust, and a cavity filled by an inviscid fluid, constituting the fluid core. This cavity is 
triaxial, which allows pressure coupling at the core-mantle boundary. The computation of 
the rotational dynamics consists in numerically integrating the Hamilton equations derived 
from the following Hamiltonian: 

'H{p,P,r,R,^,,r],,^2,V2)='Hi{P,^i,Vi,^2,V2)+'H2{p,P,r,R,^^,r],), (2) 

where "Hi is the kinetic energy of the system, and 7^2 the perturbing potential of Jupiter. 



We have (see e.g. iNovelled (120121 )) 
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and 



+ 2i/(F-iI±^)(R,- 



^)(.h*-«,6)) 



+ 2(1^ J (*^ - «f - "?) («? - "?) + i ("^^ - & - "9 («=' - "3 






+ 2W(P-^)fe-^)(OT.-66) 



5(1^ ( J (4^ - «? - •)?) ('(? - «?) + (2 - j) J (4P. - & - •,¥) {4 - &) 



+ 2.V(i'-^)(p.-^)(„,, + 6&) . (3) 



«2(p,P,r,fl,6,»h) = -|^h(:c' + !/') + f2(i' -»')), (4) 

where Q is the gravitational constant, Mq^ the mass of Jupiter, n the orbital frequency of 
lo, and d the distance lo- Jupiter, x, y and z are the coordinates of the vector pointing to 
Jupiter in the reference frame (/i, /2, /s) linked to the principal axes of inertia of lo. 

The canonical variables of the problem are: 
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P = 
R = 

Vi = 
V2 -- 


G 

--P{l-cosK), 


6 = 


= -v/2P(l-cosJ)sin/, 
= a/2Pc(1 + cos Jc) sin /c, 


= y/2P{l-COsJ)cOsl, 

= a/2Pc(1 + cos Jc) cos/c- 



(5) 



Each of these 4 hnes is related to a dynamical degree of freedom. The first one is related 
to the longitudinal motion of the whole body, p is very close to the spin angle and G is the 
angular momentum of lo. K is the obliquity with respect to the normal to the reference plane 
(here the equatorial plane of Jupiter at J2000.0), and r is the node associated. The third 
degree of freedom is related to the polar motion (or wobble) of lo, J being its amplitude. 
And the last one is related to the orientation of the velocity field of the fluid filling the core, 
we have in particular P^ = Gc/{nC) where Gc is the angular momentum of the pseudo-core. 
The pseudo-core is very close to the core, it lacks of physical relevance but is convenient to 
write a Hamiltonian formulation of the equations. These canonical variables do not directly 
represent observables of the rotation (in fact only the surface can be observed), but exact 
observable quantities can be extracted from them. 



The coordinates of Jupiter come from real ephemerides of lo, here LI. 2 (jLainey et al. 



20061 ). This way, we consider the orbital dynamics of lo with the most possible accuracy. 
The orbital period of lo is ~ 1.769 day and its eccentricity 4 x 10~^. It experiences a 462-d 
periodic perturbation due to the proximity of the 2 : 1 orbital resonance with Europa, and 
is locked into a laplacian orbital resonance involving also Europa and Ganymede. 

The numerical integrations are performed with the 10th order Adams-Bashforth-Moulton 
predictor-corrector integrat or. Once the solutions of the system have been obtained, we use 



Laskar's NAFF algorithm (lLaskarlll993l . 120051 ) to represent them as sums of quasiperiodic 



N 
n=0 

for complex variables, or 



N 



(t)^^A:exp«t) (6) 



n=0 

for real ones. In doing this, we can identify, in each variable of the problem, the influence of 
every single perturber (Jupiter, the other satellites, the Sun. . . ). 



lo is assumed to be in a dynamical equilibrium known as Cassini State 1 (jCassinilll693 
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Colombdll966[ ). as a consequence the initial conditions of the numerical integrations should 
be appropriately chosen. Deriving the equilibrium related to a simplified system (e.g. one- 
dimensional rigid rotation and circular orbit of the perturber) is usually possible, but in a 
sophisticate system as we used, with complete ephemerides, it cannot be done accurately 
enough without using a perturbation theory. It is possible to derive an approximate equilib- 
rium, but the numerical solutions will exhibit some free librations around the equilibrium, 
that are supposed to have been damped in the real system. For th is reason, we chose t o 



improve iteratively the initial conditions in using the algorithm by iNoyelles et al.l ( l2013l ) , 
consisting in 

1. starting from "pretty acceptable" initial conditions, 

2. running a numerical integration, 

3. identifying the free librations around the equilibrium, 

4. removing them from the initial condions, and reiterate the process. 

Once we have computed the equilibrum solution for our 6 models, we derive the following 
observable outputs: 

• longitudinal hbrations of the mantle 4>m = Pm — nt, often called physical librations, 
where Pm is the spin angle of the mantle, 

• obliquity of the mantle e^ (angle between the angular momentum of the mantle and 
the normal to the orbit of lo) , 

• amplitude of the polar motion J^, 

• tilt of the velocity field of the fiuid Jc, 



the formulae giving these quantities being extensively derived in lNoyelles et al.l (120 lOf ). The 
results are gathered in Tab|3l 

We can unfortunately see a high degeneracy in the sense that observing the rotation of 
lo should not allow to draw conclusions on its interior, the differences between the outputs 
being too small. The only number changing significantly is the tilt of the velocity field of 
the fluid Jc, that cannot be directly observed. 



Table 3: Variations of the outputs in the different models. 0^ is the amplitude of the 
longitudinal librations, Jm the polar motion of the mantle, Jc the tilt of the fluid, and em is 
the obliquity of the mantle. 



Models 


0m 

(462 d) 


0m 

(1.76 d) 


<Jm> 


< Jc> 


<e^ > 


Cm 

(274 d) 




(arcsec) 


(arcsec) 


(mas) 


(arcsec) 


(arcsec) 


(arcsec) 


1 


39.82462 


30.73077 


136.592 


5.89884 


7.88443 


2.35306 


3 


39.82475 


30.60912 


135.689 


5.52936 


7.87640 


2.33113 


5 


39.82475 


30.58332 


138.148 


5.41307 


7.87537 


2.32659 


2 


39.82303 


32.35491 


139.912 


15.82812 


8.14896 


3.60967 


4 


39.82346 


31.89343 


140.550 


13.13468 


8.06166 


3.08910 


6 


39.82345 


31.88809 


136.018 


12.78276 


8.04716 


3.07908 



4. Elliptical instability 



In all the calculations, we have assumed that the flow of the fluid is laminar. In fact, as 



i nitial ly seen experimentally by iPierrehumbertI (119861 ) and theoretically explained by iBayly 
( 119861 ) in the context of an unbounced strained uniform vortex, the periodic forcing of the 
elliptical cavity on the underlying rotation state produces a pairwise resonance of inertial 
waves which can grow exponentially. These st udies follow ind ependent predictions in the 
1970s, a recent review of the topic is given in iKerswelll ( 12002| ). To check the stability of 
the flow we need to cons ider the growth rate a for an arbitrary perturbation -y, see e.g. 
Kerswell fc MalkusI ( Il998l ): 



a{t) 



ld\a<v^/2> 
2 dt 



< V -VU ■v> 



where the flow U is the velocity field of the fiuid a priori assumed to be laminar, and 
<>= / d V, V being the volum e of the fluid. The flow is stable when a < 0, and unstable 
otherwise. ICebron et al.l (12012[ ) have recently derived the following formula: 



a 



n 



17 R 

— ep 
64 



2.62(1-7/)- 



T] 



16 



. 5-" ../' (9) 

in the context of a triaxial body perturbed by a primary whose rotation axis is close to 
the geometrical polar axis, and is under the influence of a magnetic fleld. The parameters 
involved in this formula are: 



e: amplitude of the physical librations (Tabj3]) at the orbital frequency (period: 1.76 
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d) 

/3 = 






rj: ratio between the internal and the external radii of the fluid layer. It is equal to 
for fully liquid cores. 

Ekman number E = ^^ where z/ is the fluid kinematic viscosity, and Q = nPc/6 = n 
the mean velocity of the fluid. This is a kind of adimensional viscosity. 



Elsasser number A = -^-^, where Bq is the intensity of the magnetic fleld and cxe the 
flu id electrical conduct ivity. This contribution o f the magnetic field ha s been derived 
by ICebron et al.l ( 12012| ) , generalizing a result by iHerreman et al.l (|2009[ ) . 



As iHerreman et al.l fl2009[ ) and ICebron et al.l f l2012l ) did, we took Bq = 1850nT, u 



lo- 



rn 



^s ^ and (Je = 4 X 10^ S.m ^. The results are given in TabJD 



Table 4: Grow rate a, and the time associated compared with the spin-up time t. 



spin— up- 



Models 


a 


1/a 


^spin—up 




{yr-') 


(kyr) 


(kyr) 


1 


1.142 X 10^3 


0.875 


3.321 


3 


1.561 X 10-=^ 


0.640 


3.140 


5 


1.700 X 10^3 


0.588 


3.109 


2 


3.325 X 10-^ 


3.008 


4.770 


4 


2.260 X 10"^ 


4.425 


4.559 


6 


2.353 X 10^^ 


4.250 


4.557 



The mean growth rate a is positive in every model. Anyway, it is smaller in the "even" 
models, i.e. with a FeS core, with a growth time of the order of 3, 000 years, while it is 
always smaller than 1,000 years in the "odd" models ( Fe core). For comparison we give in 
the last column the spin-up time tspin-up = ^/{n\/E) (JGreenspan fc Howardlll963l ). This 
is the typical spin-up/spin-down time necessary for the fluid to recover the mantle velocity. 
tspin-up is usually assumed to be long enough so that the velocity of the fluid can be considered 
as constant. We can notice that it is of the same order of magnitude as the growth time l/cx. 

In all these mo dels we flnd that th e inert ial waves of the core of lo should be un- 
stable. The work of iKerswell fc MalkusI (119981 ) showed that the waves are unstable, and 



Herreman et al.l (120091 ) conflrmed the result even when the effect of the ma gnetic fleld is 



included. This last reference suggests a growth time of ~ 63 years. However, ICebron et al. 
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( 120121 ) argue that this growth rate has been calculated in an extreme optimal case, in par- 
ticular in considering an extremum of the instantaneous departure from the synchronous 
rotation. They conclude that the inertial waves in lo's fluid core should be stable. Our re- 
sult lies between these two opposite conclusions, i.e. a positive but quite small growth rate. 
The main reason why our result is different from Cebron's with the same formula is that 
our (3 (equatorial elli pticity of the core) is bip, ger. They also assumed a smaller amplitude of 
libration, derived by 



Comstock fc Bills! ( 120031 ) for a rigid lo. 



5. Energy budget 

lo is known for its energy dissipation resulting in volcanos at its surface. Several at- 
tempts have been made to quantify this dissipation. One way to proceed is to use observa- 
tions of Jo's surface: 



Veeder et al.l (119941 ) estimated the global heat flow to be bigger than 2.bW/'w?, 
Matson et al.l (120011 ) estimated it to be smaller than ISW/m"^, 



Rathbun et al.l (12004 ) estimated to be between 2.0 and 2.6W/m'^, by studying several 



hot spots. 



the reader can find additional references in IVeeder et al.l ( 120121 ) . Tab. 2. Another way is to 
try to detect the influence of tidal dissipation on the orbit of lo b y comparing astrometric 
observations with dynamical models. Recently, iLainey et al.l (120091 1 estimated the tidal dis- 
sipation to be E = (9.33 ± 1.84) x 10^^ W, yielding = 2.24 ± 0.45 W/m'^ assuming that 
energy is transported out of lo at the same rate. 



The energy dissipation due to the instability of the fluid core is ( iLe Bars et al.ll2010f ): 



E 



6TT 



-—p,RtV^\n' 



QO|5/2 I ,5/2 

" \^so\ 



(10) 



where M^ is the mass of the fluid core. Re its radius, pc its density, f2'^ = nP the spin rate of 
the mantle, and uso is a normalized frequency of the spin-over mode, appearing in case of 
instability. We have uso = ^^ /^^ = 1/-P, since the spin rate of the fluid is assumed to be 
n. 

Using formula ( jTOj) and \Q^ — Q°\ = 2e\ cosnt\ where e ~ 0.004 is the orbital eccentricity 
of lo, we get a mean E of 8.47 x 10^ W for Model 1 and 2.32 x 10^ W for Model 2. This 
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ls consistent with the conclusions of iLe Bars et al.l (120 lOf ) announcing a peak dissipation of 
^ 4 X 10^ W. This quantity is very small with respect to the energy dissipated by tides, i.e. 
(9.33 ±1.84) X 10^3^. 

Another possible source of dissipation is radioacti ve energy, becau s e of t he decay of 
unstable isotopes in the rock component of the body. iHussmann et al.l (120 lOl ) estimate it 
between 3.08 x 10^^ W and 5.14 x 10^^ W depending on the composition of lo. 



Table 5: Energy budget of lo, at the present time. 
Source of energy 



E 



Reference 



Surface heat flow 
Tidal dissipation 
Radioactive energy 
Fluid instability 



8.33 - 10.83 X 10^3 W 

(9.33 ±1.84) X 10^3^ 

3.08- 5.14 X 10" W 

8.47 X 10^ - 2.32 x 10^ W 



Rathbun et al. (20041 

Lainev et al. (20091 

Hussmann et al. (2010) 

this note 



The results are gathered in TabJSl We can see that if the flow of the fluid constituting the 
core of lo is un stable, the ene r gy inv olved is negligible with respect to the tidal dissipation, as 
already stated iLe Bars et al.l (120101 ). So, this cannot be responsible for Jo's volcanic activity. 



Conclusion 



In this study we used the Poincare- Hough model to predict the rotation of lo, depending 
on its internal structure. For that, we elaborated 6 interior models, in considering pure iron 
or FeS compositions of the core, and the presence or not of a different crust. This rotational 
model describes the behavior of a 2-layer body composed of a rigid mantle and a fully liquid 
core, in considering pressure coupling at the core-mantle boundary. 

For each case considered, the amplitude of the short longitudinal librations is about 
30 arcsec, and the mean obliquity ~ 8 arcsec, these two quantities being a little bigger for 
a eutectic FeS core, which is larger. But the main difference is in the tilt of the angular 
momentum of the fluid constituting the core, that is bigger (between 10 and 15 arcsec) for 
the FeS than for the pure Fe one (~ 5 arcsec). These differences should unfortunately not 
be detectable. A study of the elliptical instability indicates that the inertial waves should 
be unstable in any case, the growth time being between 500 and 800 yr for a Fe core and 
between 3 and 5 kyr for a FeS core. This last calculation only considers the influence of the 
longitudinal librations, without involving the tilt of the fluid. This instability does not have 
a significant impact on the energy dissipated at the surface. 
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